SCIENTIFIC 

REPpRTS 



X 




SUBJECT AREAS: 
BIOPHYSICS 
GENETIC MODELS 
SYSTEMS BIOLOGY 
MODELLING AND THEORY 

Received 
2 March 2012 

Accepted 
11 May 201 2 

Published 
25 May 2012 



Correspondence and 
requests for materials 

should be addressed to 
F.L. (fliu@nju.edu.cn) 

orW.W. (WangWei© 
nju.edu.cn) 



Dynamic Mechanism for the 
Transcription Apparatus Orchestrating 
Reliable Responses to Activators 



Yaolai Wang, Feng Liu & Wei Wang 



National Laboratory of Solid State Microstructures and Department of Physics, Nanjing University, Nanjing 2 1 0093, China. 

The transcription apparatus (TA) is a huge molecular machine. It detects the time-varying concentrations of 
transcriptional activators and initiates mRNA transcripts at appropriate rates. Based on the general 
structural organizations of the TA, we propose how the TA dynamically orchestrates transcriptional 
responses. The activators rapidly cycle in and out of a clamp-like space temporarily formed between the 
enhancer and the Mediator, with the concentration of activators encoded as their temporal occupancy rate 
(^tor) within the space. The entry of activators into this space induces allostery in the Mediator, resulting in 
a facilitated circumstance for transcriptional reinitiation. The reinitiation rate is much larger than the 
cycling rate of activators, thereby Rtor guiding the amount of transcripts. Based on this mechanism, 
stochastic simulations can qualitatively reproduce and interpret multiple features of gene expression, e.g., 
transcriptional bursting is not mere noise as traditionally believed, but rather the basis of reliable 
transcriptional responses. 

During cellular signaling, upstream signals may be transmitted into the nucleus, inducing variations in the 
nuclear concentrations of transcriptional activators. The transcription apparatus (TA), which mainly 
consists of general transcription factors (GTFs, including TFII-B, -D, -E, -F, and -H), RNA polymerase 
II (Pol II), the Mediator, and transcriptional activators, detects such variations and accordingly initiates mRNA 
transcripts at appropriate rates 1 " 3 . This process, namely the transcriptional response that exploits the genetic 
information to modulate cellular functions, is of critical significance to life. There has been a great progress in 
exploring the mechanism of eukaryotic transcription. The primary structures of the core components of the TA 
such as Pol II and the Mediator were determined, and the basic structural organizations among different 
components of the TA were also worked out 4 " 6 . On the other hand, the TA has evolutionarily made use of the 
stochastic nature of molecular interactions to perform functions. In vivo imaging technologies such as fluor- 
escence recovery after photobleaching (FRAP) provided new insights into this dynamic process 7 " 11 . Nevertheless, 
it is still a great challenge to unravel how the TA dynamically governs transcriptional responses 7 ' 912 " 14 . 

Despite gene specificity, the general features of activator- regulated transcription can be summarized as fol- 
lows 4-612 . Transcription begins with activators binding to the enhancer (proximal or distal to the core pro- 
moter) 412 , through which specific enzymes are recruited to alter chromatin architecture 15 . GTFs and Pol II 
then bind to the matured core promoter, forming the preinitiation complex (PIC). Subsequently, the DNA 
template strand is positioned into the active center cleft of Pol II, forming the open complex (OPC). After 
preparations such as the phosphorylation of its carboxy-terminal domain, Pol II gets away into elongation 
and begins transcribing a full length of mRNA. The complex left behind at the core promoter, termed the scaffold 
complex (SCF), facilitates transcriptional reinitiation 412 . The bridge between the enhancer-bound activators and 
the basal transcriptional machinery (composed of GTFs and Pol II at the core promoter) is the Mediator, which is 
evolutionarily conserved and crucial to almost all Pol II transcription across eukaryotes 61617 . 

To probe the fundamental dynamic mechanism by which the TA orchestrates transcriptional responses, 
we begin with a minimal model focusing on the essential properties of the TA as mentioned above. This model 
consists of the core promoter, a proximal enhancer containing one binding site for the activators (distant enhancers 
can be considered as proximal after DNA looping through the tethering elements 918 ), and related proteins (Fig. la). 
Using the theories of probability and statistics, we propose the dynamic mechanism for transcriptional initia- 
tion. To further validate this mechanism, we construct a stochastic model of gene expression, and the simulation 
results qualitatively coincide with seven experimental gene-expression profiles. We find that stochastic molecular 
interactions can be maximally exploited to ensure reliable transcriptional responses to activators. 
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Figure 1 | The minimal model of eukaryotic transcriptional initiation and its properties, (a) The model. Upstream signals induce variations in the 
concentration of transcriptional activators that modulate transcriptional initiation through the Mediator. Other proteins such as histones are not 
depicted in the diagram, (b) The relationships between the transcriptional events, (c) The activators' temporal occupancy rate R™ 0R is a monotonically 
increasing function of the number of nuclear activators, n a (n a is incorporated in a on la 0 ff). Each dot represents an independent trial with m= 10 (for more 
cases, see Fig. SI). The black line denotes the average value, (d) The standard deviation of R™ 0R , SD™, vs. a on /a 0 jf for different m. SD™ decreases with 
increasing m. 



Results 

Mathematical characterization of the dynamics of the TA. The 

dynamics of the TA are dictated by how its components are 
spatially and temporally organized on the promoter. Because the 
TA can assume many distinct configuration states and the state 
evolution is essentially stochastic, involving numerous molecules 
and complex interactions, we employ the theories of statistics and 
probability to investigate the dynamics of the TA. For simplicity, we 
assume that the concentrations of transcription-associated species 
such as GTFs remain constant around the model gene and the 
molecular interactions involving the promoter are in dynamic 
equilibrium. The term "dynamic equilibrium" does not mean that 
molecular interactions are all reversible; rather, it just requires that 
the TA should retrieve its current state after some time. A model gene 
and all the species around it constitute a system. The above 
assumptions imply that such a system is in a steady state. Let us 
consider a statistical ensemble consisting of a large number of such 
essentially identical systems, with each evolving independently. 
The number of the systems is large enough so that all possible 
configuration states of the TA can be covered by this ensemble. 
That is, each state undergone by an individual gene maps the states 
of other genes in the ensemble, and the proportion of genes in a 
special state X (e.g., genes with their enhancers bound by 
activators), P(X), remains constant over time. Equivalently, if an 
individual gene is observed at any time, the probability that the 
gene is in the state X is also P(X). In this sense, which state an 
individual gene is located in is a random event. 

For the minimal model (Fig. la), we define all configuration states 
of the TA as a universal set Q and the various states with the same key 
features as the following sub-sets, respectively (Fig. lb). A denotes 
that the enhancer is bound by an activator. S denotes that the core 
promoter is bound by the proteins of the SCF. M denotes that a 
nascent mRNA is in gestation (including the process from PIC 
formation to Pol II's escape into elongation). / denotes that the 
enhancer-bound activator is conjectured to the SCF, PIC or OPC 
through the Mediator. Because eukaryotic transcriptional initiation 
requires the presence of the SCF on the core promoter 412 , MCS. 
According to the definitions, /CAS. In the set M/, M and A are 
concurrent, i.e., the enhancer-bound activators can directly affect 
Pol II's action through the Mediator. Thus, the transcriptional ini- 



tiation under direct regulation of activators is described by the set M/, 
whereas the basal, activator-independent transcriptional initiation is 
included in the set M—J. The probability of a nascent mRNA in 
gestation, i.e., the probability that an mRNA is generated, is 



P(M) = P (Aj |S) P(S)P (j| AjS) P(M| J) + q, 



(1) 



where q is a constant representing the basal transcriptional initiation, 
and Aj is a subset of A (see SI of Supplementary Information for 
details). In Aj, the enhancer-bound activators are obligate for con- 
tacting the SCF-, PIC-, or OPC-joined Mediator. Equation (1) char- 
acterizes the relationship between mRNA production and the 
dynamic properties of the TA. 

Encoding the concentration of transcriptional activators. Owing 
to distinct architectures of promoter chromatin in different tran- 
scriptional stages, the enhancer-bound activators may perform 
various functions, such as promoting histone acetylation and 
recruiting GTFs 4 ' 515 . Specifically, the set Aj involves the enhancer- 
bound activators that are responsible for handling the basal 
transcriptional machinery and controlling transcriptional initiation. 
Moreover, the activities of these activators are also associated with 
the encoding of the nuclear concentration of activators, for P(Aj|S) 
is the only factor in equation (1) depending on the concentration of 
activators. Here, we investigate the dynamics of such activators. 

Activators move rapidly within the nucleus, and the probability of 
them reaching the enhancer is proportional to their nuclear abundance 9 . 
Let us consider a time period, during which the activators involved in 
the set Aj bind to and then depart from the enhancer for m (m= 1, 2, 3, 
...) cycles. We define the temporal occupancy rate Rtor of those acti- 

vators as R™ 0R = (^=i *£») / (E/=i EJLi < # )> where tL 
and denote the binding and unbinding time of the j-th cycle, 

respectively. For the fixed number n a of activators in the nucleus, we 
have 
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where a on and a 0 jare the propensity functions of binding and unbind- 
ing, respectively (see S2 of Supplementary Information for details). a on is 
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a function of n a , whereas a 0 ffis independent of n a . Equation (2) indicates 
that as m rises, R™ 0R converges to a deterministic value, which is a 
monotonically increasing function of n a (Fig. lc-d and Fig. SI; this is 
a general property and can be applied to cases where the number of 
cognate binding sites on the enhancer is greater than one (see equations 
S13-S18)). This convergence implies that even the time-varying con- 
centration of activators can be encoded by Rtor> provided that the 
activators cycle on and off the enhancer frequently enough over a time 
window with their concentration nearly unchanged. Indeed, there are 
active disassociation mechanisms that guarantee the rapid cycling of 
activators 919 " 22 . The binding time was estimated to be within the range 
of seconds to tens of seconds 9,10 . Moreover, it was proved on the endo- 
genous CUP1 gene that such rapid cycling is functional 10 . Presumably, 
the R TO r encodes the concentration of transcriptional activators. On the 
other hand, during the time period where the activators cycle on and off 
the enhancer for m times, the probability that the enhancer is found 
bound by such an activator is R™ 0R . Because the average of R™ 0R over 
the ensemble is also j{n a ), we have 

P(Aj\S)=f(n a ). (3) 

The constraint conditions ensuring reliable transcriptional 
responses. Given the stochasticity in the occurrence of tran- 
scriptional events, to achieve a reliable transcriptional response 
requires that the Rtor code, which timely represents the 
concentration of activators, should be transduced with high fidelity 
into the amount of transcripts. Ideally, if P(S), P(j\AjS) and P(M\J) 
all equaled 1, the exact information transduction would be 
accomplished. In the following, we present the conditions under 
which those three factors can be large enough to ensure reliable 
transcriptional responses in the presence of random fluctuations 
(Figs. S2 and S3 provide intuitive explanations for this subsection). 

Equation (3) implies that the concentration of activators cannot be 
sufficiently encoded without the persistence of the SCF on the pro- 
moter. Thus, the SCF should assemble rapidly when the chromatin 
architecture allows and be much more stable than the enhancer- 
bound activators (Condition I). Such stability of the SCF was 
observed experimentally, and the binding time of TBP (TATA- 
binding protein, the core component of the SCF) on the promoter 
can be up to 20 min in human cells 11 . For P(j\AjS) , Aj is a precondi- 
tion of the occurrence of /. Because Rtor is determined by the indi- 
vidual short binding times of activators, / should happen 
immediately after the occurrence of Aj (Fig. S3). Otherwise, the 
information about Rtor is largely lost or even falsely utilized to direct 
transcriptional initiation (note that / is a precondition of M). 
Therefore, to correctly transfer the Rtor code, the Mediator should 
act by waiting to bind the cycling activators and transmit the 
information through allostery 23 25 (Condition II). This is because 
other kinds of molecular interactions like free collisions cannot pre- 
cisely convey the information about the binding time of activators. 
Such allosterism of the Mediator is supported by the previous work 26 . 
P(M\J) determines how the information about Rtor inherited by / is 
converted to guide the amount of transcripts. Because Rtor depends 
on the intermittent binding of activators, a large P(M\J) requires that 
during the short binding periods, transcripts should be produced at a 
rather rapid rate (Fig. S2) (Condition III). This feature is also verified 
by computational estimates of the experimental data (see S3 of 
Supplementary Information). Therefore, all the three conditions 
can be satisfied naturally. 

The dynamic mechanism of activator-regulated transcription. 

The above three constraint conditions together determine how the 
TA operates. There repetitively arises a state in which a relatively 
stable clamp-like space is formed between the Mediator and the 
enhancer (Fig. 2; according to Conditions I and II). As the SCF 




Figure 2 | Illustration of the dynamic mechanism for the TA 
orchestrating a reliable transcriptional response. There repetitively arises 
a state where a relatively stable clamp -like space is formed between the 
Mediator and the enhancer. Transcriptional activators cycle in and out of 
this space rapidly. Only when this space is occupied by activators do Pol lis 
initiate/reinitiate transcription at a faster rate than the cycling rate of the 
activators. 



was shown experimentally to be not very stable 11 , this space is 
temporarily constructed. The clamp-like space attracts free 
activators and then rapidly peels them off, with Rtor decided by 
the activators' concentration (according to equations (2-3)). Once 
one activator molecule gets into this space, there arises allostery in 
the Mediator, resulting in a facilitated circumstance for the GTFs and 
other related proteins to perform their functions. Consequently, Pol 
lis can initiate/reinitiate transcription very rapidly (according to 
Conditions II and III), with the Rtor governing the quantity of 
transcripts. 

This mechanism suggests that the molecular interactions invol- 
ving the promoter obey elegant dynamic principles as follows. 
Whereas the clamp -like space is temporarily formed, it is much more 
stable than the activators settled in it. The activators can cycle in and 
out of the space for many times even during short episodes when 
their concentration nearly remains unchanged; thus, the concentra- 
tion of activators can be represented by Rtor timely. Because the 
Mediator transmits the information via allostery and the rate of 
transcriptional reinitiation is much larger than the cycling rate of 
activators, the R T or c °de is effectively employed to direct mRNA 
synthesis. In a word, the clamp-like space is the structural basis for 
reliable transcriptional responses. Rather than being an obstacle, the 
stochastic nature of the molecular interactions is fully utilized to 
induce transcription reliably; this largely depends on different 
extents of the stability of the components of the TA, which span 
several orders of magnitude. The above arguments are supported 
by experimental data, and the typical time scales are as follows: the 
half lifetime of the clamp -like space is about 5 min 11 , the occupancy 
time of activators in the space is within the range of seconds to tens of 
seconds 10 , allostery usually occurs within the time of milliseconds to 
no more than 1 second 23 " 25 , and it just takes several seconds to rein- 
itiate a transcript (see S3 of Supplementary Information). 

Validation of the mechanism by numerical simulations. To further 
verify the proposed dynamic mechanism, we build a simplified 
stochastic model of gene transcription with physiologically realistic 
parameters (see Fig. S4 and S4 in Supplementary Information for 
details). This model depicts the key state transitions of the TA and also 
simply describes the related chromatin dynamics, thereby capable of 
charactering the transcriptional response to transcriptional activators. 
In the following, the "input" and "output" denote the nuclear 
concentration of activators and the amount of gene products, 
mRNA or protein, respectively. 
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Figure 3 | Transcriptional responses to activators based on the proposed dynamic mechanism. The input equals a on la 0 ft which is positively related to 
the nuclear concentration of activators, (a) Temporal evolution of the number of cellular mRNAs in a single diploid cell with different input levels. 
mRNAs produced by two alleles are shown separately in red and black. The transcriptional burst becomes dense with increasing the input strength, 
(b) The average input/output relationship in individual diploid cells. The maximal outputs are normalized to 1. Error bars denote the standard deviation 
of the output, SD 0Ut . The inset shows the ratio of SD out to the mean output vs. the input. Because the abundance of mRNAs or proteins also depends on 
their degradation/inactivation rates, which are modulated by cellular signaling, the rate of mRNA production more directly reflects the dynamics of the 
TA (see also Fig. S9, where the production rate of proteins is also shown 44 ' 45 ), (c) The curves of the SD 0Ut vs. the input. These curves nearly remain bell- 
shaped even at various degradation rates of mRNAs or proteins (see also Fig. S10). (d) The distribution of mRNA levels across a cell population for 
different input levels. The bin size is 10. (e) The state evolution of a promoter in response to a periodically varying input. G x denotes the enhancer is bound 
by an activator. SCF denotes the core promoter is bound by the SCF. OPC denotes the core promoter is in the OPC state. The curves describe the input, the 
corresponding states of the promoter, and the production of mRNAs (from top to bottom), respectively, (f) ChIP simulation of the transcriptional 
response. The input and the symbols are the same as in panel (e). TAT An and Pol II denote that the core promoter is bound by histones and Pol II, 
respectively. 



First, we explore the temporal evolution of the number of cellular 
mRNAs at constant input levels (Fig. 3a). Notably, mRNAs are pro- 
duced in a burst-like manner, consistent with the prevailing view 14,27 " 32 . 
For low-level inputs, one allele is transcribed while the other is silent in a 
diploid cell, and thus the bursting phenomenon is apparent. For high- 
level inputs, however, both the two alleles burst frequently so that the 
sum becomes almost constant. This suggests that the phenotype of 
persistent elevated transcriptional responses may be observed at high 
input levels 14 . 

A recent experimental analysis precluded the possibility that the 
chromatin environment plays a central role in shaping transcrip- 
tional bursting 32 . Here, we demonstrate that a burst of transcripts 
originates from persistent reinitiation by Pol lis when the clamp-like 
space is occupied by the activators (Fig. 4). That is, the initiation of 
mRNA is itself burst-like. The bursting is not mere noise; instead, it is 
a direct manifestation of the Rtor code, which represents the con- 
centration of activators and guides mRNA production. 

Second, we investigate the average input/output relationship of 
transcriptional response. The average output resembles a Hill func- 
tion of the input, which is widely used in systems biology to model 
gene expression 3,33,34 (Fig. 3b). The curve of the standard deviation 
SD 0Ut of the output versus the input is approximately bell-shaped 
(Fig. 3c). The intensity of the intrinsic noise, defined as the ratio of 
SD out to the mean output 35 , is inversely correlated with the input (the 
inset of Fig. 3b). Moreover, the above features are insensitive to the 
slight fluctuations in the input (i.e., extrinsic noise) (Fig. S5), suggest- 
ing the robustness of the transcriptional response to noise. All these 
results are in good agreement with the experimental measurements 
in both Saccharomyces cerevisiae 36 and Drosophila embryos 37 . 
Particularly, the left side of the SD 0Ut curve is lower than its right 
side; this feature is almost quantitatively in accordance with 



the experimental data 37 (see S5 of Supplementary Information for 
further discussions). By contrast, deviations from the dynamic prin- 
ciples proposed above (including the circumstances under which the 
activators' cycling is slow, the scaffold complex/clamp -like space is 
not stable, or/and the rate of transcriptional reinitiation is low) would 
reduce the capability of the TA to respond reliably to the input (Fig. 
S6). 

The input/output relationship observed in Drosophila embryos 
was believed to be realized by maximally utilizing the limit of 
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Figure 4 | The essence of transcriptional burst. Shown is a microscopic 
view of a transcriptional burst. 'CA' denotes that an activator is in the 
clamp-like space. 'OPC denotes that the transcription machinery is in the 
OPC state (a zoom-in panel is also displayed). When an activator molecule 
is present in the clamp-like space, rapid transcriptional reinitiation results 
in a burst of mRNAs. 
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molecular interactions 37 " 39 . The properties of such microscopic inter- 
actions are integrated to be macroscopically manifested as SD 0Ut . The 
SD out curve is still bell-shaped overall compared with the SD in curve 
(cf. Fig. Id). That is, the signature of the Rtor c °de can be directly 
transmitted to the output. This confirms that the temporal occu- 
pancy rate of activators is really exploited to regulate transcription 
and the Mediator transmits the information through allostery. On 
the other hand, the SD 0Ut curve is asymmetrical, with the right side 
being higher than the left. The reason is obvious. When the input is 
very high, the enhancer is bound by the activators almost all the time, 
and thus the fluctuations mainly reflect the dynamic properties of the 
SCF and the transcriptional reinitiation by Pol lis. Our further simu- 
lations show that the right side of the SD 0Ut curve drops as the 
stability of the SCF or the rate of transcriptional reinitiation is 
increased; only when increasing their strength beyond the physio- 
logical ranges can the curve become symmetrical (Fig. S6F). This also 
verifies that both P(S) and P(M\J) are indeed large enough. 
Therefore, the properties of SD 0Ut should conclusively prove the 
microscopic transcriptional mechanism. 

Third, we probe the distribution of mRNA levels across a large 
cell population exposed to the same input (Fig. 3d). For small inputs, 
the bursting phenomenon is especially obvious, and most cells have 
no or few mRNAs. This is consistent with the experimental obser- 
vation 27 ' 30,31 . But the distribution gradually becomes normal as 
the input increases. For a on /a G ff >1, the distribution becomes 
sharper with increasing the input. These results await experimental 
identification. 

Fourth, we simulate the transcriptional response to a periodically 
varying input. The microscopic process on a promoter is rather 
dynamic and stochastic, with different components of the TA exhib- 
iting distinct stabilities (Fig. 3e). However, the amount of mRNAs 
can follow the input. These results are in good agreement with the 
results revealed by FRAP, i.e., the TA is a highly dynamic appar- 
atus 8101 1>22 . On the other hand, chromatin immunoprecipitation 
assays (ChIP) simulations, which characterize the temporal evolu- 
tion of the distribution of different states of the promoter among a 
cell population, reveal a strong regularity in the distributions 
(Fig. 3f). The patterns of both the activators binding to the enhancer 
and the SCF and Pol II binding to the promoter follow the input. 
mRNA transcripts are produced in phase with the input, whereas 
histones occupy the promoter in a reversed phase. All these results 
are well consistent with the experimental findings 22,40 . Therefore, the 
observed discrepancies between the results from FRAP and ChIP 
experiments may originate from different resolutions involved in 
the measurements. ChIP measurements integrate molecular interac- 
tions both temporally and over the cell population, whereas FRAP 
more tightly reflects the instantaneous interactions. Moreover, the 
transcriptional response to the time -varying input is robust to 
extrinsic noise but sensitive to composite input signals, and devia- 
tions from the dynamic principles (such as the cases with low cycling 
rate of activators, the unstable scaffold complex/clamp -like space, 
or/and low rate of transcriptional reinitiation) would reduce the 
response capability (see Figs. S7 and S8). 

Discussion 

Based on the general features of the eukaryotic TA, we have proposed 
the fundamental dynamic mechanism by which the TA orchestrates 
reliable transcriptional responses to cellular signals. Although our 
work is built on the general architecture of the TA, different profiles 
of gene expression can be accounted for simultaneously. This implies 
that the eukaryotic TA likely shares the same basic mechanism in 
mediating transcriptional responses, just as the same set of GTFs is 
involved. We have shown that the TA is an elegant apparatus; the 
stabilities of its components are widely differentiated such that the 
stochastic nature of molecular interactions is employed to achieve 
reliable transcriptional responses. The activators' cycling in and out 



of the clamp-like space modulates the amount of mRNA transcripts 
initiated through the Mediator's allostery. The concentration of 
activators are represented by the statistical quantity Rtor> which 
essentially leads to burst-like mRNA production. Thus, the tran- 
scriptional bursting is not only the phenotype but also the basis of 
reliable transcriptional responses. 

Traditionally, it was believed that activators regulate gene tran- 
scription through recruiting proteins such as GTFs and Pol II 41 . Here, 
we argue that this is mainly the process of TA assembly. Our results 
demonstrate that it is through controlling the circumstance where 
Pol lis reinitiate transcription that activators mediate the responses 
to upstream signals; the clamp-like space between the enhancer and 
the Mediator is the structural basis for Rtor to guide the amount of 
mRNA transcripts. It is worthy to note that the transcriptional mech- 
anism proposed here can be viewed as a complicated realistic version 
of the clustering model 13 ' 27,42 ' 43 . In that model, a promoter is in an "off 
or "on" state, and only when the promoter is in the "on" state, the gene 
is transcribed actively. Here, we further show that Pol lis reinitiate 
transcripts rapidly only when the clamp-like space is occupied by 
activators. Finally, it is worth mentioning that it is still difficult to 
experimentally elucidate the dynamic mechanisms of molecular 
machines, whereas our approach employing statistics and probabil- 
ity theory is powerful in this field. 

Methods 

Mathematical derivations. The relationships between the transcriptional events 
depend on the basic knowledge that eukaryotic transcriptional initiation requires the 
SCF 412 and only the enhancer-bound activators rather than free activators can 
regulate transcription. Equation (1) was derived using the probability theory. 
Equation (2) was derived by employing the Gillespie theory 46,47 . The implications of 
equations (1-3) and the constraint conditions were also intuitively illuminated in 
Figs. S2-S3. All mathematical derivations are detailed in Supplementary Information. 

Stochastic simulations. The stochastic model was constructed based on the proposed 
dynamic mechanism of activator-regulated eukaryotic transcription and related 
experiments. The Gillespie algorithm 47 was used to perform simulations. The very 
detailed description of the model, the parameter values, and the methods of numerical 
simulations are presented in Supplementary Information. 
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